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ABSTRACT 

We present stellar evolution models for 0.5 - 1.2 Mq at scaled metallicities of 
0.1 - 1.5 Zq and 0/Fe values of 0.44 - 2.28 O/Poq . The time dependent evolu¬ 
tion of habitable zone boundaries are calculated for each stellar evolution track 
based on stellar mass, effective temperature, and luminosity parameterizations. 

The rate of change of stellar surface quantities and the surrounding habitable 
zone position are strong functions of all three quantities explored. The range of 
orbits that remain continuously habitable, or habitable for at least 2 Gyr, are 
provided. The results show that the detailed chemical characterization of exo¬ 
planet host stars and a consideration of their evolutionary history are necessary 
to assess the likelihood that a planet found in the instantaneous habitable zone 
has had sufficient time to develop a biosphere capable of producing detectable 
biosignatures. This model grid is designed for use by the astrobiology and ex¬ 
oplanet communities to efficiently characterize the time evolution of host stars 
and their habitable zones for planetary candidates of interest. 

Subject headings: astrobiology, (stars:) evolution, catalogs, (stars:) planetary 
systems 


1. INTRODUCTION 

One of the fundamental reasons for planet searches, and indeed one of the central mo¬ 
tivations in the held of astronomy, is the eventual discovery of life on a planet outside our 
solar system. Though stellar evolution is not the only component to consider when working 
to understand planetary habitability, it is nonetheless one of the most important and phys¬ 
ically well-understood factors. It plays a huge role in creating a “habitable” environment 
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by providing energy to orbiting planets. Current astrophysical research into the habitabil¬ 
ity of exoplanets focuses mostly on the concept of the classical habitable zone (HZ), the 
range of distances fr om the star over which liquid water could exist on a planet’s surface 
flKasting et al.l 1199311 . The location of the HZ is determined on the stellar side primarily by 
the host star’s luminosity and secondarily by its spectral characteristics. These properties 
serve as boundary conditions for a planetary atmosphere calculation that predicts a planet’s 
surface temperature, and therefore the possibility for stable liquid water. 

Over time, our understanding of the HZ has become more sophisticated as models of 
planetary atmospheres and their interaction with incoming stellar radiation have advanced. 
There are several components we must consider in order to accurately gauge the surface 
conditions of a planet, including: planetary atmospheric composition, cloud cover, and 
whether certain features (e.g. H 2 O absorption bands) could be detected in atmospheric 


spectra fjvon Paris et ahl 120111) : the physics of CO 2 cloud condensation and evaporative 
water-loss; geophysical exchange with an atmosphere; pl anet mass, density, and surface 


temperature; and assumptio ns about greenhouse effects flKasting et al.l Il993l: ISelsis et ah 


2 OO 7 I: iKopparapu et ahl 1201411 . For this work we use reasonably conservative estimates for 
the inner and outer boundaries of the HZ, based on HZ limit equations derived from a 
radiative-convective, 1-D atmospheric code that includes updated m odels of water absorp¬ 
tion in planetary atmospheres, discussed in IKopparapu et ahl (120141) (see §2.3). 


It has become commonplace to announce the discovery of a planet when it is found in 
the HZ of its host star. A test of instantaneous habitability requires over-plotting a planetary 
orbit on the inner and outer boundaries of the HZ corresponding to observed temperature 
and luminosity of the star. This is not as straightforward as it would appear, since uncertain¬ 
ties in the stellar luminosit y (often fro m distance uncertainties) can cause the HZ to move 
more than its entire width (iKanel 1201411 . Nonetheless, for stars with well-measured parame¬ 
ters we can place a planet within the HZ, and we can estimate the number of planets that 
should be in the HZ in a statistical sense. The Kepler mission was launched in 2009 with an 
explicit goal of hnding an Earth-sized planet orb iting within the HZ , and was very success¬ 


ful in discovering t hese small terrestrial objects flBorucki et al.ll201ll: ICharpinet et al.ll2011 


Fressin et al.ll2012ll. Subsequent data ana . 


( Catanzarite fc Shaol 2011 : Petigura et al 


ysis has predicted that one in Eve Sun-like stars 
2OI3I: ^idosIboiO l and approximately half of all 


M-dwarf stars f Borucki et al. I2OI0I: Batalha et al. 2013[) are likely to host an Earth-sized 


planet in the HZ. 


For astronomers hoping to hnd planets with detectable biosignatures, this introduces 
the novel problem of having too many candidates. Directly detecting and obtaining spectra 
for an Earth-like planet around a Sun-like star is typically difficult and would likely require 
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a large-scale space mission. Such a n iission may be able t o observe no more than about 
a hundred stars over its lifetime (e.g. iTurnbull et ahl (120121) h As a hrst step, it would be 
useful to narrow-down the planetary candidates and associated host stars on the basis of 
whether a planet is currently habitable, which we can estimate directly from the measured 
properties of the star; however, the current habitability of a planet is not the only important 
consideration for hnding Earth-like worlds. If the ultimate goal is the discovery of life, the 
evolution of stars and the associated HZ must be taken into account as well. 


The likelihood for detection of life on a habitable planet ultimately stems from a prob¬ 
ability distribution constructed of multiple terms falling into three broad classes: the ability 
of a planet to sustain a biosphere, the likelihood that biosignatures can persist at detectable 
levels on the planetary surface, and the technique being used to observe the planet. We can 
begin dehning individual terms for each of these categories, which will be expanded upon 
as our understanding progresses. As examples of the latter classes, some inhabited planets 
will probably be unidentihable until we are able to examine them from up close. An obvious 
example is Europa in our own solar system, which may support a biota under its icy crust. 
The limited communication of the potentially life-sustaining region with the surface gives 
Europa a low potential for detectability because biosignatures will escape at low levels at 
best and will be rapidly destroyed by Jupiter’s radiation environment. 


More subtle effects may come into play that could even prevent detection of surface bio¬ 
spheres on terrestrial planets. As examples, planets with a highly reducing bulk chemistry 
could conceivably £x free oxygen faster than even a substantial population of photosynthe- 
sizing organisms could produce it, which could deprive us of an opportunity for detection 
of life. Free oxygen is often suggested as one of the more “useful” or “indicative” biosig¬ 
natures, although recently that idea has been called into question; it has become more 
apparent that oxygen p resent in a planetary atinosph e re may not necessarily indicate any 
sort of biogenic proces s (jPomagal-Goldman et al.ll2014J: IWordsworth fc Pierrehumbertll2014J: 




huger fc BarnesI I 2 OI 5 I) . Such a planet might, on the other hand, be able to sustain high 


levels of biogenic methane that could not be attributed to geological processes such as ser- 
pentinization. This planet would have a higher likelihood of detection for missions that are 
able to observe methane as well as oxygen, which is why a thorough characterization of the 
entire planetary system is important. 


One of the hrst things we should consider in attempting to understand a whether a planet 
may support life, and the likelihood of detection, must be the characteristics of the host star 
and the location of the surrounding HZ. We have the soundest basis for developing detection 
strategies for Earth-like life, and such life is clearly favored by Earth-like conditions. Such 
conditions are, by dehnition, most likely to be found in the HZ. However, the instantaneous 
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location of the planet and HZ are too limited to assess this aspect of the planet’s habitability. 
The location of the HZ is a function of time; we know that as stars age, the habitable region 
will move away from the star due to the gradual increase in stellar luminosity (L). The 
effective temperature {Teff) will also evolve, changing the spectrum impinging upon the 
planet’s atmosphere. Therefore, any orbiting planets detected around a particular star may 
only have spent a relatively short time in the HZ. It is easy to assume that a planet around 
a sedate, hve billion year old star will be as rich and diverse as Earth. Time dependent 
models show, however, that more than half of the orbits that are in the HZ at some point 
during the evolution are actually only habitable during the latter part of the star’s life. Our 
hypothetical hve billion year old planet will not be a very fecund place if it just entered 
the HZ 100 Myr ago. It is estimated that life only produc ed a detectable change in Earth’s 
atmospheric chemistry 1-2 Gyr after Earth’s formati on flKastindll993l: iBrocks et al.lll999 


Kopp et al.ll2005l: lAnbar et al.l 120071: 1 Crowe et al.ll2013h even though Earth was in the Sun’s 
HZ from very early times. 

Identifying these kinds of systems from readily measurable stellar and orbital parameters 
would rapidly reduce the number of habitable planetary candidates. An essential parameter 
that affects stellar and HZ evolution is the detailed chemical composition of a star. The 
rate of stellar evolution and the change in Teff an d L are dependent on the abundances 
of individual elements, especially oxygen and iron flYoung et al.N2012l) . The sensitivity of 
stellar evolution t o composition arises from two effects: the eq uation of state (EOS) and the 
radiation opacity flRoeers et al.lll996l: llglesias fc RogersI 119961) . The changes in the EOS are 
relatively minor, but rearranging the proportions of diffe rent species at a consta nt [Fe/H] 
can result in signihcant opacity changes of tens of percent (llglesias fc Rogerslll9961) . Opacity 
changes the rate of leakage of radiation, and increased radiation pressure in the stellar 
envelope drives expansion, resulting in larger radii and lower effective temperatures. A slower 
rate of energy loss also requires a slower rate of nuclear burning to maintain hydrostatic 
equilibrium. So we expect stars with enhanced abundance ratios [X/Fe] to be cooler, less 
luminous, and longer-lived relative to stars with the same [Fe/H]. 

Measurements of [Fe/H] alone are insufficient to predict the stellar and HZ evolution 
flYoung et al.N2012h . In practice, it is usually only the iron abundance that is measured 
for many stars, and the rest of the elements are assumed to scale in the same proportions 
present in the Sun; this means models of stellar evolution for stars of different metallicity 
are generally created under this assumption. This is the nearly universal practice in stellar 
modeling, though the abundance ratios in real stars vary substantially. Except for a uniform 
enhancement of the a elements (O, Ne, Mg, Si, S, Ar, Ca, Ti) in stars of very low metallicities, 
variations in abundance ratios at a given [Fe/H] are neglected, despite the fact that stellar 
evolution is sensitive to specihc composition. 
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For solar mass stars near solar [Fe/H], in all cases, enhanced abnndance ratios cause 
larger changes than depleted ratios. C, Na, and Mg have small but noticeable effects. C has 
a high abundance but relatively few electron transitions and a low ionization potential. Mg 
and Na have lower abundances but higher opacity per gram than C, resulting in a similar 
degree of shift in the tracks. Si has less impact due to its smaller range of variation, and Al, 
Ca, and Ti have very small effects due to their small abundances; these can be neglected for 


the p urpose of habitability. The largest changes arise from variation in oxygen (lYonng et a,l. 


20121 ). Both L and T^fj of the enriched compositions are systematically lower at a given age, 
but stellar evolution models show the most profound effect is on the pace of the evolution. 
For example, for a 1 Mq star with solar [Fe/H] the main sequence (MS) turnoff for a model 
representing the low end of the distribution of oxygen abundance in nearby stars occurs at 
an age of ~ 9 Gyr. Solar composition has a turn-off age of ~ 10 Gyr, and an oxygen-rich 
model turns off at ~ 11 Gyr (and at lower L and T^ff). 

In this paper, we examine how stellar mass, metallicity, and individual elemental abun¬ 
dance ratios influence the long term habitability of planetary systems. We present a grid 
of stellar evolution models for Sun-like stars, with masses of 0.5 - 1.2 Mq , metallicities 
of Z = 0.1 - 1.5 Zq , and a spread in oxygen values ranging from 0/Fe = 0.44 - 2.28 
O/Fe© . When calculating the time-dependent HZ inner and outer radii for each of our 
stellar evol utionary tracks, we hay e chosen to use a widely cited and well-known set of pre¬ 


scriptions flKopparapu et al.ll2014|) . though the main utility of our catalog is that our models 


are easily transferrable for use with any particular set of HZ limit equations; this includes 
scenarios where a so-called “cold start” is not viable (i.e. a frozen planet enters the HZ 


KoDDaranu et al. 

.7 

(2014 

Following 

Kastine et al. 


work under the assumption that a cold start is in fact plausible. 
(119931) we present a scenario where the outer limit of the HZ at the 


zero age main sequence (ZAMS) is a hard limit that does not vary with the star’s evolution 
so that these two scenarios can be compared. 


We examine F, G, and K-type stars in particular because they are the closest in physical 
properties to our own Sun, which is a good starting point of comparison if we want to look 
for potential host stars of habitable Earth-like planets. M-stars present other issues for 
habitability, including high levels of stellar activity with high energy particle and x-ray 
fluxes, close-in HZs with potential for tidal locking, spin-orbit resonances and tidal heating, 
and dynamically packed inner systems. We therefore defer calculations of M-stars, which 
require more comprehensive treatments of complex molecular opacities, to a separate paper. 


Ultimately, we would like to learn what kind of star could potentially host a planet that 
has remained in a continuously habitable zone (GHZ) for at least 2 Gyr. This is approxi- 
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mately the amount of time it took for life on Earth to change the atmospheric composition 
sufficiently that there would be a detectable biosignature if viewed from another star system, 
utilizing missions recommended in the most recent Decadal Review of Astronomy and Astro¬ 
physics (e.g. transmission spectroscopy with JWST or direct detection with a coronagraph, 
starshade, or interferometer). We dehne a CHZ for the star’s entire MS lifetime, as well as 
the range of distances that would be continuously habitable for at least 2 Gyr at some period 
in the star’s evolution. We describe our choice of parameter space, stellar evolution code, 
and assumptions for calculating the HZ in §2, our interpretation of the results of the models 
in §3, and our conclusions in §4. 


2. METHODS 
2.1. Parameter Space 

In this work we present a grid of stellar models suitable for the prediction of HZ loca¬ 
tions. The most important variable for stellar evolution is, of course, mass. Second is total 
metallicity, Z. Third is the oxygen abundance ratio, 0/Fe. The ratios C/Fe and Mg/Fe also 
produce a small ef f ect in stellar evolution and have substantially variable a bundance ratios 
( Neves et al. 2009 : Mishenina et ah 2008 : Takeda 2007 : Young et ah 2014). These will be 
explored in a future paper. In this work, ratios without brackets (e.g. 0/Fe) indicate the 
linear absolute abundance ratio in terms of mass fraction, while a bracketed ratio denotes 
the log of the atom number relative to the solar abundance value for that same element. 
The latter is the conventional [0/Fe] given by 


logio 


jO/Fe) 
iO/Fe] 


o 


= |0/Fe] 


( 1 ) 


This work primarily quotes linear ratios relative to solar (i.e. 0/Fe = 2.28 O/FeQ ) since 
the range of abundance ratios is small enough to not require logarithmic notation. We use 
mass fraction as this is the conventional usage for stellar evolution calculations. 

Here we consider the major contributors, mass, metallicity, and oxygen abundance. 
Variations in Z alone are made with a hxed abundance pattern that is uniformly scaled. 


Gonzalez Hernandez et al. 2010: Hinkel et al. 


2007: 

Bond et al. 

2006, 

2008 


20141) . The values of O/Poq are taken from 


Young et al.l (120141) . which analyzes the intrinsic spread (not accounted for by observational 
error) in elemental abundance ratios [X/Fe] for 5 large surveys. We use the [0 /Fe] from the 


median survey, which is consistent with the values used for the calculations in lYoung et al. 
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(120121) ■ A solar composition from iLoddersI (120101) was adjusted to the mean abundance ratios 
of the median sample for the elements observed, with all other elements being maintained 
at solar values for standard composition. The adjusted values are listed in Tabled] Changes 
in 0/FeQ at each metallicity are made by changing the absolute abundance of O while 
holding all other metal abundances constant. The abundances of H and He are adjusted in 
compensation to ensure the sum of mass fractions = 1. 


Determinations of O abundance are notoriously sensitive to non-LTE effects (iGrevesse et al 


20071 ) and line-blending, which may introduce random unphysical errors in the O abundance 
measurement that are not accou nted for in the quoted observational errors. The surveys 
considered in I Young et ahl (120141) were chosen in part for using the most reliable of the op¬ 
tical O lines, but we also examine a smaller range o f O variation. Instead of choosing a 
random value, we use the Uintrinsic (lYoung et al.ll20141) for a sample of 40 “solar twins” from 
Ramirez et al.l (12009[) . Qualifying as solar twins requires observable parameters quantita¬ 
tively close to the Sun. In order for a star to match the Sun closely in physical observables, 
it is necessary for it to have a composition close to solar. Using Uintrinsic for the solar twin 
sample thus provides a conservative estimate of [0/Fe] variation. 


The initial grid encompasses solar-type stars on the MS, covering a mass range from 
0.5 - 1.2 Mq at increments of 0.1 Mq , corresponding to spectral types MO - FO at Z = Zq . 
Models are calculated for metallicities of 0.1 - 1.5 Zq at increments of 0.1 Zq . Two models 
are calculated with values of O/Feg at each Z value (end members 2.28 and 0.44 0/FeQ ). 
The other two O cases (0.67 and 1.48 0/FeQ ) were done only at solar Z value. The grid is 
complete for the MS until hydrogen exhaustion in the core for all cases. Post-MS evolution 
for short lived stars (> 1 Mq ), cool stars (e.g. M-dwarf stars), and contributions from minor 
elemental constituents (such as C and Mg) will be explored in a future publication. 


2.2. TYCHO 


The stars included in our catalog were simulated using the stellar evolution code TYCHO 
(lYoung fc ArnettI 120051) . TYCHO is a 1-D stellar evolution code wi th a hydrodynamic for¬ 


mulation of the stellar evolut i on equations. It uses OPA L opacities (jlglesias fc Rogerslll996 


Alexander fc Ferguson! 119941: iRogers fc NavfonovI l2002h. a combined OPAL and T immes 
equation of state (HELMHOL TZ) (Timmes fc Arnett 1999 ; Rogers fc Navfonov 2002 ). grav¬ 
itational settling (diffusion) (IThoull et al.l Il9941) . general relativistic gravity, time lapse, 
curvature, automatic rezoning, and an adaptable nuclear reaction network with a sparse 
solver. A 177 element network terminating at ^^Ge is used throug;hout t he evolution. 
The network uses the latest REACLIB rates (iRauscher fc ThielemaunI l2000l: I Angulo et al. 

























































1999 

Iliadis et al. 

2001 ; 

Wiescher et al. 

2006) 


fennnfi 


, wea k rates from lT;anga.nke fc Martmez-Pinedo 
), and screening from iGraboske et al.l fllQTsI) . Neutrino cooling from plasma processes 
and the Urea process is included. Mass loss is include d but is trivial for this mass range. 


(Heightened early mass loss seen in some young stars flWood et al.l 120051) is not included.) 


It incorporates a description of tur bulent convection f Meakin fc Arnett 


2007: Arnett et ah 


20091. 12 OIOI: lArnett fc Meakinll201111 . which is based on 3-D, well-resolved simulations of con¬ 


vection between stable layers, which were analyzed in detail using a Reynolds decomposition 
into average and fluctuating quantities. It has no free convective parameters to adjust, unlike 
mixing-length theory. 


TYCHO outputs information on stellar surface quantities (including L and Tgjf) for 
each time-step of a star’s evolution, which we then use to calculate the inner and outer 
radii of the HZ as a function of the star’s age. The initial composition of the stellar models 
were adjusted as described in §2.1, and new OPAL opacity tables were generated for the 
specihe abundances for each 0/Fe value to match the composition of the stellar model. 
TYCHO begins calculation with a fully convective model on the Hayashi track. We limit 
our discussion of the HZ to the span between the Zero Age Main Sequence (ZAMS) and 
the Terminal Age Main Sequence (TAMS). For this purpose, the ZAMS is dehned as the 
luminosity minimum coinciding with the beginning of complete hydrogen burning in the core. 
The TAMS is dehned as the time of hydrogen exhaustion in the core {Xh < 1.0 x 10“®). 
The complete tracks are publicly available o 


2.3. Calculating Habitable Zone Extents 


TYCHO evolutionary tracks are used to estimate the extent of the 
i n the stella r evolution. For thes e estimates we follo w th e prescriptions of 


dZ at each point 


Kopparapu et al. 


fj2013l. l2014(l , which follow from ISelsis et al.l fl2007|) and iKasting et al.l fjl993l) . These pre¬ 
scriptions parameterize the orbital radii of the HZ as a function of L and T^ff facilitating 
translation from stellar evolution tracks to HZ distance estimations. 


Kopparapu et ahl (1201411 use radiative-convective planetary atmosphere models with in 


put synthetic stellar spectra produced by the PHOENIX code (lAllard fc Hauschildtl 119951) . 
The predicted distance from a star for both the inner and outer edges of the HZ is param¬ 
eterized as a function of stellar temperature and luminosity, for several kinds of planetary 
atmospheres. The calculations do not include other sources of heating such as tides (small 
for these stars) or effects such as eccentric orbits. The distance d of the HZ boundary is 


^http://bahamut.sese.asu.edu/~payoung/AST_522/Evolutionary _Tracks_Database.html 
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defined by 


d 


L/LeY^ 

Seff ) 


( 2 ) 


where Seff is an effective flux received by the planet based on the radiative transfer 
calculations for different initial spectra corresponding to different Teff passing through the 
planetary atmosphere. 


Seff = Seffe + aT, + bT^ + cT^ + dT^ (3) 

Where T* = Teff - 5780 K, SeffQ is the effective flux for Earth from the Sun at the solar Teff, 
and a, b, c, and d are coefhcients to a polynomial £t. Each distinct planetary atmosphere 
model generates a unique set of coefficients. 


Five cases are described in iKopparapu et al.l (120131) : (1) Recent Venus, (2) Runaway 
Greenhouse, (3) Moist G reenhouse, (4) Maximum Greenhouse, and (5) Early Mars. However, 


Kopparapn et al.l (120141 ) replaces the Moist Greenhouse case entirely with the Runaway 


Greenhouse case, since the differences between the two are minimal. We us e the updated 
values for the coefficients (a, b, c, d) and Se//© from iKopparapn et al.l (120141 ). 


The four remaining cases represent a signihcant range of HZ approximations when we 
consider both the stellar flux incident on the planet due to the planet’s distance from the 
host star, along with several different atmospheric properties that may exist on an Earth- 
like planet. Though the amount of solar radiation the planet receives is the largest factor in 
determining HZ distances, the wavelength-dependent radiative transfer and radiative losses 
from atmospheres with different levels of greenhouse gases, water, and clouds must also be 
considered. HZ boundaries are undoubtedly strongly influenced by the presence of clouds in 
a planet atmosphere. We would expect H^O cloud s to move the inner HZ boundary inward 
(lKastindll988l: ISelsis et al.l 120071: lYang et al.l 120131) because their contribution to a planet’s 
albedo generally outweighs the contribution to the greenhouse effect. Gonversely, GO2 ice 
clouds are expected to cause warming in a dense GO2 atmosphere because they reflect thermal 
radiation back to the planet’s surface more e fhciently than they reflect incoming radiation 
back to space (IForget fc PierrehumbertI 119971) . 


The Runaway Greenhouse (RGH) and Maximum Greenhouse (MaxGH) cases represent 
the conservative estimates for the inner and outer radii of the HZ, respectively. The RGH 
limit is the distance at which a planet’s oceans would evaporate entirely. The inner edge 
of the HZ is determined by the level of water saturation in the planet’s atmosphere (a 
warm environment causes evaporation of H2O from the planet’s surface) and the subsequent 
rapid loss of hydrogen to space due to heating effects from the proximity of the host star. 
Gonversely, planets that orbit near the outer HZ boundary could develop dense, G02-rich 
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atmospheres through outgassing. CO 2 begins to condense out of the atmosphere at a certain 
distance away from the parent star (due to colder temperatures) which reduces the overall 
greenhouse effect. CO 2 is also an effective Rayleigh scatterer (2.5 times better than air), and 
so a dense CO 2 at mosphere is ex pected to have a high albedo which would further offset the 
greenhouse effect flKastindllOOll ). Thus, the conservative outer HZ boundary (MaxGH) is 
the location where Rayleigh scattering by CO 2 begins to outweigh the greenhouse effect. 

Additionally, the Recent Venus (RV) case predicts an inner HZ edge much closer to the 
star, while the Early Mars (EM) case predicts a more distant outer boundary. These we 
will refer to as the optimistic cases. The RV case limit from radar observations of Venus 
by the Magellan spacecr aft that suggest liquid water has been absent from the surface of 
Venus for at least 1 Gyr fISolomon fc HeadlllQQll) . when the Sun was about 90% as luminous 
as it is today. This gives a more optimistic empirical estimate for the inner HZ boundary. 
Likewise, the (optimistic) Early Mars (EM) outer boundary has been estimated based on 
the observation tha t early Mars was wa rm enough for liquid water to flow on its surface 


( Pollack et ahl 119871: Bibring et a 


raised some debate fjSegura et ah 


2006). Al though the issue of a warmer early Mars has 


2 OO 2 I. I 2 OO 8 I) . this gives a good estimate of a more distant 


outer boundary at which a planet could potentially remain habitable. 

The TYCHO evolutionary tracks are used as input to CHAD (Calculating HAbitable 
Distances), which is a code we dev eloped to calculate the inner and outer HZ boundaries for 


each of the HZ limits presented in iKopparapu et ahl (120141) . CHAD is easily upgradable to 


incorporate improved HZ predictions as they become available in the future. In this paper 
we will focus on the conservative estimates, but all of the calculations will be made publicly 
available in machine-readable format along with the evolutionary tracks. 


RESULTS 


From a stellar point of view, three main factors influence the time evolution of the HZ. 
These are L and Te//, their rate of change, and the stellar MS lifetime. We conhrm th a t mas s 
and total metallicity influence these factors considerably. Following on lYonng et ahl (120121) . 
the ratio of 0/Fe is also signihcant over the entire range of M and Z. When we consider the 
evolution of the HZ, it is clear that the instantaneous HZ calculated from observed stellar 
properties is often not a good indicator of a high likelihood for detection of biosignatures, 
except for very low mass stars. Detectability ultimately depends on assumptions made about 
the timescales involved in evolution of life and biogeochemical evolution of biosignatures. 
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3.1. Stellar Properties and Main Seqnence Lifetimes 

TYCHO stellar evolution tracks are used to determine the stellar parameters of interest 
at each time-step in a star’s evolution. The evolution was calculated from the pre-MS to 
H exhaustion. The pre-MS evolution is rapid, with a similar timescale to planet formation, 
and can be neglected. This provides us with a MS lifetime and a rate of change for the 
stellar properties. The Hertzsprung-Russell diagrams (HRD) shown in Figures [T] and [2] 
demonstrate the effect of variations in oxygen abundance ratios for each mass in our data 
set, from 0.5 M© to 0.8 Mq in Figured! and from 0.9 Mq to 1.2 Mq in Figure [2l For both 
hgures, the solid lines represent solar-value for the 0/Fe ratio, whereas the dashed lines 
and dotted lines correspond to depleted and enriched oxygen cases (0.44 O/FeQ and 2.28 
O/Fe© ), respectively. Signihcant changes are seen when the O abundance ratio varies, even 
at constant Z. Trends with 0/Fe parallel those with the total metallicity; low 0/Fe stars are 
bluer, more luminous, and shorter-lived. Higher 0/Fe stars are cooler, less luminous, and 
longer lived. Thus, the specihc O abundance in stars (not just the overall metallicity) plays 
a signihcant role in stellar evolution and, consequently, planetary habitability. 

Figure [3] shows log(L/L 0 ) vs. time (Gyr) for the compositional end member cases in our 
library for a 0.5 Mq star. Models for 0.1 Zq , 1.5 Zq , Zq , 0.44 O/Fe© and 2.28 O/Fe© (each 
at Z = Zq ) are shown. Figures 0] and [5] show the same for 1 Mq and 1.2 Mq stars. For 
all masses we see a signihcant variation in the MS lifetime, with the lowest metallicity star 
living just two thirds as long as the highest metallicity model for the 0.5 Mq star, and about 
half as long as the highest metallicity model for the 1 Mq star. Of particular importance 
is the change wrought by varying 0/Fe. In this case, an otherwise solar composition star 
with 0/Fe = 2.28 0/FeQ actually lives longer than a Z = 1.5 Zq star. This illustrates the 
danger in only measuring [Fe/H] and assuming that a scaled solar metallicity will accurately 
predict a star’s evolution. For the 0.5 Mq star, the high-end cases for O and Z are quite 
similar in terms of L, however the enriched O case sees a longer MS lifetime by about 6 Gyr. 
As expected, we see that stars with higher Z (or 0/Fe) are less luminous, but considerably 
longer-lived. Table 0] shows MS lifetime estimates for the range encompassed by our catalog, 
acquired using this method. 

A more subtle effect to consider is the rate of change of Te// and L. Low opacity models 
undergo a larger change in luminosity than do the higher opacity models at the same mass, 
over a shorter total lifetime. Thus, dL/dt is greater for low metallicity (or 0/Fe) models, 
especially during the second half of a star’s evolution on the MS. The radial movement of the 
HZ boundaries is concomitantly faster. At 0.5 Mq , a 1.5 Zq model increases in L/Lq by 
0.025 dex, while a 0.1 Zq model is brighter by 0.05 dex at the end of the MS. This translates 
to the luminosity of the star changing by a factor of 5.4 at 0.1 Zq compared to only 2.9 at 1.5 
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Zq . The absolute change in luminosity is even more sensitive to variations in composition at 
higher masses, though the percentage change is smaller. Tables |3] and H] show the total change 
in L in terms of L/LzamSi and the total change of T^ff in units of Kelvin (respectively) 
from ZAMS to TAMS. 

Finally, because of how the slopes of dL/dt (and to a lesser extent dT/dt) change over 
time, the range of orbits in the HZ at different points in the MS evolution can change 
substantially. When considering the potential for detectability, we wish to avoid planets 
that have only recently entered their HZ. Table 0 shows the fraction (listed as percentages) 
of orbital radii which enter the HZ after the midpoint of the MS. It is clear from these results 
that a third to a half of orbits that are in the HZ only become habitable in the second half 
of the star’s life. The effect is more pronounced at higher mass and composition. 


3.2. Location of the Habitable Zone 

Figures [6] and [7] show the HZ from a perspective perpendicular to a hypothetical orbital 
plane, demonstrating how the HZ can vary between different stars. Figure [6] shows a 0.5 
Mq star and a 1.2 Mq star (end member masses) at the lowest and highest metallicity cases 
(0.1 Zq and 1.5 Zq , respectively). HZ boundaries are solid for the ZAMS and dashed for the 
TAMS. The high mass stars exhibit some overlap between the outer edge at the ZAMS and 
the inner edge at the TAMS, which might correspond to a “Continuously Habitable Zone” 
(see §3.3). The largest overlap appears for the 0.1 Zq model. Given this HZ prescription, 
there are no orbits around the low mass stars that remain within the HZ for the entire MS. 
Figure [7] similarly shows the HZ distance evolution from ZAMS to TAMS for the extreme 
oxygen cases (2.28 0/FeQ and 0.44 0/FeQ ) at Zq . We see that the HZ can change sub¬ 
stantially over the MS depending on the host star’s detailed composition. Tables |6] and [7] 
show changes in the location of the HZ radius in AAU from ZAMS to TAMS for both the 
inner boundary (RGH) and the outer boundary (MaxGH), respectively. 

Although there is nearly as much variation in the location of the HZ for low and high 
mass stars and low and high metallicity, the rate of that change has very different implications 
for habitability. As demonstrated in Table [21 the least massive star’s MS lifetime is much 
longer (ranging from ~ 60 - 99 Gyr) than any reasonable timescale for the development of a 
detectable biosphere. Most orbits even near the boundaries will remain in the HZ for billions 
of years. Conversely, a more massive star’s entire MS lifetime at low metallicity is only ~ 
2.8 Gyr. A planet that takes 2 Gyr or so to evolve a detectable biosphere (like Earth) would 
need to remain in the HZ for nearly the entire MS of the host star. Given the rate of change 
in the HZ position, very few orbits would qualify. 
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Figure [8] shows the inner and outer edges of the HZ for each stellar mass for all com¬ 
positions at the ZAMS, while Figure [9] shows the same information at the TAMS. Higher 
opacities result in boundaries at lower radii. It is interesting to note that with increasing 
mass, there seems to be a widening of the overall HZ range, as well as a larger spread in 
the HZ distances due to compositional variation. Figures |8] and |9] also show the spread in 
O variation for hve values (0.44, 0.67, 1, 1.48, and 2.28 0/FeQ ) at solar metallicity only. 
These models are indicated by elongated solid lines. The scaled Z cases and the varied O 
cases appear to be consistent in that they both show the HZ spreading out with increasing 
mass. The range in distance of the HZ edges for the oxygen values at Zq is smaller than 
for the entire range of metallicity. This is expected, since the total change in opacity of 
the stellar material is much larger for a factor of hfteen change in total metallicity than a 
factor of two change in oxygen abundance. Note that these hgures also include the spread 
in oxygen calculated at each metallicity value. The elongated dotted lines represent the end 
member values for the spread in oxygen abundance (0.44 and 2.28 O/Fe© ) calculated at end 
member metallicity values (0.1 and 1.5 Zq ). These models extend the range of HZ distance 
even further than do the models for oxygen abundances at Zq alone. The difference in HZ 
position as a function of composition is larger for higher mass stars because the absolute 
change in L is larger for higher masses. The position of the outer edge changes more than 
that of the inner because the behavior of the Maximum Greenhouse scenario is more sensitive 
to the spectrum of the incoming radiation and therefore Tg//. 

We produce complete evolutionary tracks for the position of the HZ as a function of 
time for all stellar models. With an independent age estimate for the star, as well as its 
mass and its composition, the position of an extrasolar planet can be compared not only to 
the current HZ, but also its past and future location. Assuming stellar properties are well 
measured, the dwell time of an observed exoplanet in the HZ can be estimated to the level of 
accuracy of the atmosphere models predicting HZ boundaries. Figure [TUI shows the location 
of the inner (solid) and outer (dashed) edges of the HZ as a function of time for three stellar 
masses (0.5, 1, and 1.2 Mq ) at Zq and high, standard, and low 0/Fe values. The smallest 
radii correspond to the most enhanced oxygen model. 

Again, Table [2] provides MS lifetimes for each stellar model of interest. The least massive 
(0.5 Mq ) star has the same pattern of MS lifetimes for each of the three oxygen models 
that we see with the higher mass stars, though the MS lifetimes are much longer for the 
0.5 Mq star. The enriched oxygen model (2.28 0/FeQ ) is estimated to hve 99.3 Gyr, while 
the depleted oxygen model (0.44 0/FeQ ) is estimated to hve 83.5 Gyr, making the spread 
in lifetimes about 15.8 Gyr. Gonversely, the highest mass star (1.2 Mq ) has the shortest 
overall MS lifetime for each of the oxygen models, with a spread of only about 1.2 Gyr 
between the two end-member oxygen cases. However, because the average MS lifetime of 
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the more massive star is shorter, the percent difference in MS lifetime estimates between the 
end-member oxygen cases is much larger for the 1.2 Mq star than for the 0.5 Mq star. 

Each of the stellar models represented in Table [2] generally have longer lifetimes as Z 
increases, except for the 0.5 and 0.6 Mq stars, which show a turnover in MS lifetimes between 
1 Zq and 1.5 Zq (at enriched oxygen). The percent differences between the two values is 
4.85% for the 0.5 Mq star and 1.25% for the 0.6 Mq star. The reason for the turnover in 
lifetimes for the lowest mass stars at the highest metallicity/oxygen values is because other 
effects of high heavy element content become large enough to compete with enhanced opacity. 
In these cases, the amount of hydrogen in the core is reduced by several percent while the 
fraction of helium has increased. The increase in He reduces the Thompson scattering opacity 
of the inner, mostly ionized, region of the star while simultaneously reducing the number of 
free particles contributing to pressure support. Coupled with the reduced amount of fuel, 
the MS lifetime ends up being smaller. At even higher metallicity this trend would extend 
to higher masses. These factors working in tandem produce the turnover we observe, in that 
the MS lifetime is actually shorter for the 1.5 Zq case than it is for the 1 Zq case, both at 
enriched oxygen. 

With more oxygen present in the host star - that is, the host star would have a higher 
0/Fe ratio than present in the Sun - the HZ will be closer in to the star because the star is 
less luminous and is at a lower effective temperature, and the stellar lifetime will signihcantly 
increase. Likewise, a lower oxygen abundance will produce shorter overall MS lifetimes with 
HZ distances that are markedly farther away from the host star. The total MS lifetime varies 
by about 3 Gyr between the end-member oxygen abundances. The rate of change of L and 
Te// is much faster for a star with a shorter MS lifetime, and therefore the location of the HZ 
changes much more quickly. The habitable lifetime for a terrestrial planet varies by about 
4.5 Gyr. Figure [TUI also demonstrates that the average lifetime for a solar mass star at solar 
composition (both metallicity and oxygen value) is about 10 Gyr, as expected. 

Figure [TT] shows a similar trend for the end-member metallicity cases, for the same 
three masses (0.5, 1, and 1.2 Mq ). We hnd here that Earth would be interior to the HZ at 
all times for the lowest metallicity value. The variation is somewhat less between the high 
metallicity case and standard Z value, compared to the high-end oxygen and standard O 
value; however, they are much larger between the low metallicity case and standard value 
compared to standard and low oxygen. This is unsurprising, considering that a factor of 
2.28 increase in O is large enough to produce an opacity increase of a similar magnitude to 
an overall 150% uniform scaling of metallicity. Reducing O to 0.44 O/Fcq is a much smaller 
reduction than an overall 90% reduction in all opacity producing elements. 

A key consideration is how changing oxygen alone compares to scaling the entire metal- 
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licity of a star. For both types of compositional variation, we see the same trends in HZ 
distance with enhanced or depleted compositions. As overall scaled Z increases, MS lifetime 
increases while L decreases due to greater opacity within the star, which reduces the rate at 
which energy can escape. Similarly, as the O abundance increases, MS lifetime increases, L 
decreases, and the HZ will be located nearer to the host star. This is the same overall effect 
that metallicity exhibits, and for the same reason - increased stellar opacity. What has not 
been appreciated before is that at a given [Fe/H] other individual elements can vary enough 
to signihcantly affect the stellar evolution. In fact, the higher O case at Zq prolongs the MS 
lifetime a bit more than the high scaled Z case alone. 


3.3. Continuously Habitable Zones 


It is useful to quantify the dwell time of a planet in the HZ as a function of its orbital 
radius. It is clear that the instantaneous habitability of a planet is insufficient to determine its 
likelihood of having life or its detectability. A planet must remain habitable for a considerable 
amount of time. The simplest exercise is to hnd a continuously habitable zone (CHZ), a range 
of orbital radii that remain in the HZ for the entire MS. Figure [12] shows the CHZ for stars of 
all masses in our range at solar composition. The CHZ is dehned by considering the boundary 
overlap between ZAMS and TAMS. Low mass stars have no CHZ for the conservative cases. 
For the optimistic case, we have overlap for the entire mass range, though it is much smaller 
for low mass stars than for high mass stars. This would seem to indicate that low mass stars 
would have a low likelihood of detectability, but this is somewhat misleading due to the long 
MS lifetimes of low mass stars. 


This simple version of the CHZ is of limited utility. A more useful measure would be 
the zone that is continuously habitable for enough time for a surface biosphere to produce 
a measurable and identihable change in the planet’s atmosphere. On Ear th this process 
took ~ 2 Gyr ( Summons et ah 1999 : Kasting fc Catling 2003 : Holland 2006 ). Here we show 
results for orbits that remain habitable for at least 2 Gyr (CHZ 2 ). This assumes that Earth’s 
timescale for biological modihcation of the atmosphere is representative. This is not meant to 
suggest that other timescales are impossible, but in narrowing down a large pool of candidates 
there is an advantage in looking for what we know works from observing Earth’s history. It is 
important for us to implement some kind of screening process that allows us to narrow down 
where any upcoming planet-hnding missions should look for potentially habitable planets. 
It is possible of course that life may exist on non-Earth-like planets that do not fall under 
our 2 Gyr criteria. Figure [13] shows the CHZ 2 for standard composition. This is determined 
from the inner edge of the HZ at 2 Gyr after the beginning of the MS and the outer edge 
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2 Gyr before the TAMS, which represents the location around the star that a planet could 
remain habitable for at least 2 Gyr. Because of the much longer lifetimes of the low-mass 
stars, there is a higher proportion of the HZ included in the GHZ 2 than for the basic GHZ. 
Likewise, due to the shorter lifetimes of the more massive stars, it would be statistically 
less likely to hnd a planet orbiting this type of star that has been in the GHZ for at least 2 
Gyr. We would be less conhdent that a planet located outside of the GHZ 2 would produce 
detectable biosignatures than one within the GHZ 2 . With the exception of planets orbiting 
M-stars (because of their extremely long lifetimes), it is a useful exercise to determine the 
location of the GHZ 2 . Table |8] shows the fraction of time a planet would spend in the GHZ 2 
vs. time it would spend in the HZ over its entire MS lifetime. This information will help 
to quantify the kinds of stars we should focus on in the continued search for habitable (and 
possibly inhabited) exoplanets. 


In our consideration of the GHZ, we must also address the issue of “cold starts.” The 
previous discussion assumes that as the HZ expands outward due to the effects of stellar 
evolution, any planets that were initially beyond the boundaries of the HZ could potentially 
become habitable as so on as the HZ reac h es the m; indeed, the albedos used in the planetary 
atmosphere models of iKopparapu et ahl (120141) are relatively low, which assumes a planet 
could become habitable upon entering the HZ. However, it may be unlikely that a completely 
frozen planet (a “hard snowball”) entering the HZ late in the host star’s MS lifetime would 
receive enough energy in the form of stellar r adiation to reverse a global glaciation, especial ly 
if the planet harbors reflective GO 2 clouds fiGaldeira fc Kastindll992l: iKasting et al.lll993l) . 


Figures [TH [TSl and [16] offer an alternative scenario to that of Figures [101 [HI and [TS] 
respectively, in which we have assumed a cold start is possible. Here we treat the outer 
boundary of the HZ at the ZAMS as a hard limit that does not co-evolve with the host star 
over time, so that only planets in the HZ from the beginning of the MS will be habitable; 
if a planet is able to dwell within the HZ from early times, then a cold start is no longer a 
problem. Figure [HI shows the evolution of the HZ for the same cases as Figure [TO] now with 
cold starts prohibited. Figure [TSl demonstrates what the GHZ might look like if we £x the 
outer boundary at the ZAMS value. In this hgure, the thin-lined shaded region represents 
the range of orbits in which a planet would simply be in its host star’s HZ from the beginning 
of the MS, while the thick-lined shaded region represents the range of orbits that would be 
continuously habitable for the star’s entire MS lifetime. For low mass stars there are no 
continuously habitable orbits, though again the long lifetimes of these stars make that less 
of an issue. In all cases, however, the range of orbits that could produce an inhabited planet 
is restricted compared to the case in which cold starts are allowed. 


Similarly, Figure [16] shows the GHZ 2 when cold starts are prohibited, where the shaded 
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region is now representative of the range of orbits that would be continuously habitable for at 
least 2 Gyr (as opposed to the entire MS lifetime), and orbits that also would be in the host 
star’s HZ from the beginning of the MS, which would thus help us rule out any planets that 
would not enter the HZ until later in the star’s MS lifetime. If we attempt to understand 
a GHZ in this way, we see that the ideal habitable regions do not include the TAMS outer 
limit, where a planet could potentially start out very cold and only enter the HZ late in its 
lifetime. Gonsidering the time dependent evolution of the HZ is even more important in this 
case, since a signihcant fraction of planets that are observed to be within the nominal HZ 
now are likely to be cold start cases. 


4. CONCLUSIONS 


Given the classical dehnition of the Habitable Zone, which assumes a surface biosphere 
that supports liquid water, the properties of the host star are obviously of fundamental im¬ 
portance. However, “habitable” does not automatically mean inhabited, which in turn does 
not equate with observable biosignatures. Many factors must be considered in calculating 
the likelihood for detection. The astrophysics should be treated in a more nuanced fashion, 
and the time evolution of the location of the HZ must also be considered. 


Time evolution of the parent star is an important consideration, because the timescale 
for life to develop to a point where it alters the planetary atmosphere sufficiently for biogenic 
non-equilibrium species to be detectable may be quite long. On Earth, this process took ~ 
2 Gyr. Some attention has been paid to this aspect, but it is overlooked surprisingly often 
when announcements are made of planets in the HZ of the associated host star. Additionally, 
the concept of a GHZ is too simple, since it does not take into account the variable lifetimes 
of stars. There is no GHZ for a solar composition 0.5 Mq star (using the conservative HZ 
boundary estimates), but a set of orbits need only be habitable for a small fraction of the total 
stellar lifetime (which is on the order of 10^^ years) to meet the 2 Gyr criterion. We propose 
a 2 Gyr continuously habitable zone as an aid to estimating the likelihood for detection, 
which is a necessary, but not completely sufficient, condition. Because of the shape of the 
luminosity vs. time curve for a star, a third to a half of habitable orbits only become so in 
the second half of the star’s life on the MS. For example, one can certainly imagine planets 
in these orbits entering the HZ when they are geologically dead, which may be unfavorable 
to life (see Table [S]). 


In turn, the evolution of a star depends on its elemental composition. This is a basic 
fact of stellar astrophysics and has been considere d in HZ evolution in terms of total metal- 
licity by some groups (e.g. iDanchi fc LopezI (120131) b However, it is important to distinguish 
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between the metallicity of a star as measured by [Fe/H] and the actual abundances of indi¬ 
vidual elements. Metallicity is often used interchangeably with [Fe/H], which expresses the 
amount of iron alone relative to the Sun’s elemental abundance. The standard approach to 
stellar modeling is typically to scale the metallicity assuming all elements scale in the same 
proportions found in the Sun, though observations show that this is very often not the case 
in reality. The HZ distance can be substantially affected even when only abundance ratios 
are changed. Evaluating habitability by modeling stellar and HZ evolution requires models 
that span a range of variation in abundance ratios, as well as total metallicity. For the same 
reason, characterizing a system requires measurements of multiple elemental abundances, 
not just [Fe/H]. 

Here we present the results of changing 0/Fe at levels similar to those seen in nearby 
stars as well as the total metallicity. Both changes influence the evolution by changing the 
stellar opacity and therefore have similar effects. For low metallicity or low 0/Fe at a given 
[Fe/H], MS lifetimes are shorter, and the total luminosity change over the MS is larger. This 
results in a high dL/dt and dTgff/dt, and a correspondingly rapid change in the location 
of the HZ. What is not widely appreciated is that changes in abundance ratios can have 
very large effects. For example, variation of 0/Fe values from a base solar composition by 
an amount observed in some planet host candidates can change the stellar lifetime more 
than increasing the metallicity by 50%. An increase in O increases the total heavy element 
abundance, but there are two important differences. First, the opacity will be different for 
different mixes of elements with the same integrated heavy element abundances. Second, 
as long as model comparisons with individual stars are made based only on the measured 
[Fe/H], which is the common practice, very large errors may be present. 

Since many targets of radial velocity planet searches have high quality spectra that can 
be used to determine abundances, it should not be difficult to compare to stellar models with 
more accurate compositions as long as such models exist. We have made the library of stellar 
evolution models and predicted HZs upon which this paper is based available to the commu¬ 
nity at http://bahamut.sese.asu.edu/~payoung/AST_522/Evolutionary_Tracks_Database.html. 
A routine for interpolation between tracks for any specihed model contained within the grid 
boundaries will also be provided. The library will be extended in the future for additional 
elements and masses, as well as later evolutionary stages for sufficiently short-lived stars. 
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Table 1. Abundance adjustments from solar composition to average sample composition. 


Element 

[X/H] standard 

C 

0.173 

0 

0.051 

Na 

0.068 

Mg 

-0.029 

A1 

-0.124 

Si 

0.078 

Ca 

0.027 

Ti 

0.072 


^Standard composition 
values relative to solar 
value dodders 2010) 


Table 2. 


MS lifetimes (Gyr) at standard and end member oxygen values, at solar and end 


member metallicity values, for all masses. 


Composition 

0.5 M0 

0.6 M0 

0.7 M0 

0.8 M0 

0.9 M0 

1.0 M0 

1.1 M0 

1.2 M0 

0.1 Zq , 0.44 O/FeQ 

60.4 

33.6 

19.8 

12.3 

8.0 

5.5 

3.9 

2.8 

0.1 Zq , O/FeQ 

64.2 

35.9 

21.2 

13.1 

8.6 

5.8 

4.1 

3.0 

0.1 Zq , 2.28 O/Fe0 

76.9 

43.3 

25.6 

15.7 

10.3 

6.9 

4.8 

3.5 

Zq , 0.44 O/FeQ 

83.5 

52.8 

32.5 

20.6 

13.5 

9.1 

6.2 

4.5 

Ze I O/Fe0 

88.4 

57.3 

35.8 

22.7 

14.8 

9.9 

6.7 

4.9 

Z0,2.28 O/Fe0 

99.3 

64.4 

41.2 

26.1 

16.9 

11.1 

7.5 

5.7 

1.5 Z0,0.44 O/Fe0 

87.6 

57.6 

36.5 

23.3 

15.3 

10.3 

7.0 

5.1 

1.5 Zq , O/FeQ 

93.6 

60.8 

39.0 

24.9 

16.3 

10.9 

7.4 

5.5 

1.5 Z0,2.28 O/Fe0 

94.6 

63.6 

41.2 

26.2 

17.0 

11.1 

7.9 

5.8 













Table 3. A{L/ Lzams) for each mass and end-member composition. 


Composition 

0.5 Mq 

0.6 Mq 

0.7 Mq 

0.8 Mq 

0.9 Mq 

1.0 Mq 

1.1 Mq 

1.2 Mq 

0.1 Zq , 0.44 O/FeQ 

5.64 

4.62 

3.70 

2.99 

2.48 

2.04 

1.69 

1.43 

0.1 Zq , O/Fe0 

5.41 

4.41 

3.51 

2.84 

2.35 

1.94 

1.60 

1.37 

0.1 Zq , 2.28 O/FeQ 

5.27 

4.40 

3.50 

2.73 

2.28 

1.87 

1.52 

1.27 

Zq , 0.44 O/Fe0 

3.27 

3.21 

2.73 

2.19 

1.79 

1.50 

1.13 

1.00 

1 O/Fbq 

3.17 

3.16 

2.76 

2.21 

1.78 

1.46 

1.09 

1.01 

Zq , 2.28 O/Fe0 

3.07 

3.00 

2.71 

2.19 

1.74 

1.39 

1.01 

1.02 

1.5 Zq , 0.44 O/FeQ 

2.96 

2.98 

2.64 

2.13 

1.75 

1.47 

1.05 

0.94 

1.5 Zq , O/Fe0 

2.89 

2.84 

2.57 

2.11 

1.71 

1.41 

1.01 

0.94 

1.5 Zq , 2.28 O/FeQ 

2.81 

2.75 

2.45 

1.99 

1.58 

1.28 

1.05 

0.91 


Table 4. AT (K) for each mass and end-member composition. 


Composition 

0.5 Mq 

0.6 Mq 

0.7 Mq 

0.8 Mq 

0.9 Mq 

1.0 Mq 

1.1 Mq 

1.2 Mq 

0.1 Zq , 0.44 O/FeQ 

719 

555 

315 

136 

58 

101 

217 

168 

0.1 Zq , O/Fe0 

718 

592 

344 

149 

53 

29 

125 

114 

0.1 Zq , 2.28 O/FeQ 

652 

665 

425 

171 

82 

22 

-5 

41 

Zq , 0.44 O/Fe0 

320 

502 

495 

303 

118 

-15 

-137 

-243 

•^O 1 O/FeQ 

303 

463 

511 

346 

170 

30 

-100 

-215 

Zq , 2.28 O/Fe0 

256 

349 

439 

346 

202 

78 

-57 

-193 

1.5 Zq , 0.44 O/FeQ 

277 

417 

475 

335 

170 

44 

-128 

-244 

1.5 Zq , O/Fe0 

240 

345 

430 

330 

193 

76 

-74 

-198 

1.5 Zq , 2.28 O/FeQ 

248 

296 

344 

269 

144 

50 

-126 

-237 


Table 5. Fraction (%) of radii which 

are in the HZ after the midpoint of the MS. 

Composition 

0.5 Mq 

0.6 Mq 

0.7 Mq 

0.8 Mq 

0.9 Mq 

1.0 Mq 

1.1 Mq 

1.2 Mq 

0.1 Zq , 0.44 O/FeQ 

30.6 

33.1 

35.4 

37.9 

40.8 

44.3 

48.2 

52.0 

0.1 Zq , O/Fe0 

31.1 

33.8 

36.4 

39.0 

41.8 

45.1 

49.0 

52.5 

0.1 Zq , 2.28 O/FeQ 

31.4 

34.5 

37.1 

40.0 

42.8 

46.1 

49.7 

53.5 

Zq , 0.44 O/Fe0 

37.7 

39.1 

42.1 

45.1 

48.1 

51.1 

55.0 

56.1 

Zq , O/FeQ 

38.1 

39.2 

42.0 

45.3 

48.5 

51.9 

56.2 

56.4 

Zq , 2.28 O/Fe0 

38.4 

39.6 

42.2 

45.8 

49.5 

53.3 

58.2 

56.6 

1.5 Zq , 0.44 O/FeQ 

39.4 

40.2 

42.9 

46.1 

48.9 

52.0 

56.6 

57.3 

1.5 Zq , O/Fe0 

39.5 

40.7 

43.1 

46.4 

49.8 

53.1 

57.7 

57.8 

1.5 Zq , 2.28 O/FeQ 

40.0 

41.1 

43.8 

47.5 

51.5 

55.4 

56.5 

58.0 
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Table 6. AAU for each mass and end-member composition at inner HZ limit (RGH). 


Composition 

0.5 M0 

0.6 M0 

0.7 M0 

0.8 M0 

0.9 M0 

1.0 M0 

1.1 M0 

1.2 M0 

0.1 Zq , 0.44 O/FeQ 

0.39 

0.49 

0.59 

0.69 

0.77 

0.82 

0.83 

0.86 

0.1 Zq , O/FeQ 

0.37 

0.46 

0.55 

0.63 

0.71 

0.77 

0.79 

0.82 

0.1 Zq , 2.28 O/Fe0 

0.34 

0.42 

0.50 

0.57 

0.65 

0.70 

0.72 

0.73 

Zq , 0.44 O/FeQ 

0.23 

0.30 

0.35 

0.40 

0.45 

0.49 

0.50 

0.56 

Zq , O/FeQ 

0.22 

0.28 

0.34 

0.38 

0.42 

0.46 

0.46 

0.54 

Zq , 2.28 O/Fbq 

0.20 

0.26 

0.31 

0.36 

0.39 

0.409 

0.403 

0.51 

1.5 Zq , 0.44 O/Fe0 

0.21 

0.27 

0.32 

0.37 

0.41 

0.45 

0.44 

0.50 

1.5 Zq , O/Fcq 

0.19 

0.25 

0.31 

0.35 

0.38 

0.41 

0.41 

0.48 

1.5 Zq , 2.28 O/Fe0 

0.18 

0.24 

0.29 

0.32 

0.35 

0.37 

0.42 

0.46 


Table 7. AAU for each mass and end-member composition at onter HZ limit (MaxGH). 


Composition 

0.5 M0 

0.6 M0 

0.7 M0 

0.8 M0 

0.9 M0 

1.0 M0 

1.1 M0 

1.2 M0 

0.1 Z0,0.44 O/Fe0 

0.68 

0.86 

1.03 

1.20 

1.34 

1.42 

1.47 

1.55 

0.1 Zq , O/Fcq 

0.65 

0.80 

0.96 

1.11 

1.24 

1.34 

1.37 

1.45 

0.1 Zq , 2.28 O/Fe0 

0.61 

0.74 

0.88 

1.00 

1.13 

1.22 

1.26 

1.28 

Zq , 0.44 O/FeQ 

0.42 

0.53 

0.63 

0.71 

0.79 

0.87 

0.88 

0.99 

Zq , O/FeQ 

0.40 

0.51 

0.60 

0.68 

0.74 

0.80 

0.81 

0.95 

Zq , 2.28 O/FeQ 

0.37 

0.47 

0.57 

0.63 

0.68 

0.72 

0.71 

0.90 

1.5 Zq , 0.44 O/Fe0 

0.38 

0.49 

0.58 

0.65 

0.72 

0.79 

0.78 

0.89 

1.5 Zq , O/Foq 

0.36 

0.46 

0.55 

0.62 

0.68 

0.73 

0.72 

0.85 

1.5 Zq , 2.28 O/Fe0 

0.34 

0.44 

0.52 

0.58 

0.63 

0.67 

0.75 

0.83 


Table 8. Fraction (%) of time spent in 

CHZa vs. 

the entire MS. 


Composition 

0.5 M0 

0.6 M0 

0.7 M0 

0.8 M0 

0.9 M0 

1.0 M0 

1.1 M0 

1.2 M0 

0.1 Z0,0.44 O/Fe0 

85.8 

81.0 

75.2 

69.0 

61.6 

54.3 

O 

b 

0.0 

0.1 Zq , O/FeQ 

86.8 

82.6 

77.5 

71.5 

64.6 

57.3 

50.0 

0.0 

0.1 Zq , 2.28 O/Fe0 

88.7 

85.1 

81.0 

75.9 

69.7 

63.0 

55.8 

0.0 

Zq , 0.44 O/Fcq 

91.5 

89.2 

86.2 

82.6 

78.2 

73.0 

67.4 

59.5 

Zq , O/FeQ 

91.8 

89.9 

87.2 

83.9 

79.9 

75.3 

70.7 

62.2 

Zq , 2.28 O/Foq 

92.7 

90.9 

88.8 

86.0 

82.5 

78.2 

74.4 

65.2 

1.5 Zq , 0.44 O/Fe0 

92.0 

90.0 

87.6 

84.5 

80.3 

75.5 

71.2 

63.4 

1.5 Zq , O/Fcq 

92.7 

90.8 

88.5 

85.6 

81.8 

77.0 

73.1 

65.1 

1.5 Zq , 2.28 O/Fe0 

92.8 

91.3 

89.4 

86.8 

83.3 

78.7 

70.6 

65.7 


'No orbits are continually habitable for 2 Gyr as a result of the short MS lifetime. 
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Fig. 1.— HRD, Evolutionary tracks from ZAMS to TAMS for 4 masses, with 0/Fe = 0.44 
O/Fe 0 (dashed), 1.0 O/Fe© (solid), and 2.28 O/Fe© (dotted) all at Z = Zq . The rightward- 
most dotted line is for the 0.5 Mq star with enriched oxygen, while the leftward-most dashed 
line is for the 0.8 Mq star with depleted oxygen. 
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Fig. 2.— HRD, Evolutionary tracks from ZAMS to TAMS for 4 masses, with 0/Fe = 0.44 
O/Fe 0 (dashed), 1.0 O/Fe© (solid), and 2.28 O/Fe© (dotted) all at Z = Zq . The rightward- 
most dotted line is for the 0.9 Mq star with enriched oxygen, while the leftward-most dashed 
line is for the 1.2 Mq star with depleted oxygen. 











24 



Fig. 3.— Log(L/L 0 ) vs. time (Gyr) for a 0.5 Mq star for five different compositions. The 
total MS lifetime varies from 65 Gyr to nearly 100 Gyr. The shortest lifetime corresponds 
to a star with total metallicity Z = 0.1 Zq . The longest lifetime corresponds not to the 
highest scaled Z model, but rather to the model with 0/Fe = 2.28 O/Fe© (at Z = Zq ). 
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Fig. 4.— Log(L/LQ ) vs. time (Gyr) for a 1 Mq star for five compositions. We see now 
that the total MS lifetime has a much smaller variation, from ~ 6 - 11.5 Gyr. Again, the 
shortest lifetime corresponds to a star with metallicity Z = 0.1 Zq , and the longest lifetime 
corresponds to the model with 0/Fe = 2.28 O/Fe© (at Z = Zq ). 
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Time (Gyr) 


Fig. 5.— \jOg{L/ Lq ) vs. time (Gyr) for a 1.2 Mq star for five compositions. Again, the 
shortest lifetime corresponds to a star with metallicity Z = 0.1 Zq , and the longest lifetime 
corresponds to the model with 0/Fe = 2.28 0/FeQ (at Z = Zq ). 
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Fig. 6.— Habitable Zone Ranges: 0.5 Mq (top), 1.2 Mq (bottom), for 0.1 Zq (left), 1.5 
Zq (right). HZ shown at ZAMS (solid) and TAMS (dashed). Inner boundaries represent 
Runaway Greenhouse and outer boundaries represent Maximum Greenhouse. 
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Fig. 7.— Habitable Zone Ranges: 0.5 Mq (top), 1.2 Mq (bottom), for 0.44 O/Fe© (left), 2.28 
O/Fe 0 (right), at Zq . HZ shown at ZAMS (solid) and TAMS (dashed). Inner bonndaries 
represent Rnnaway Greenhouse and outer boundaries represent Maximum Greenhouse. 
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Habitable Zone Distance (AU) 

Fig. 8 .— Inner and outer HZ edges (RGH and MaxGH) for all masses and compositions 
at ZAMS. The elongated solid lines represent the 5 oxygen cases at Zq for each mass: 2 
enriched oxygen models (1.48 and 2.28 0/FeQ ), 2 depleted oxygen models (0.67 and 0.44 
O/Fe© ), and 1 model representing the standard oxygen case (solar 0/Fe). The elongated 
dotted lines represent the end member values for the spread in oxygen abundance (0.44 and 
2.28 O/Fe 0 ) calculated at end member Z values (0.1 and 1.5 Zq ). From left to right are 
the highest to lowest opacity models. It is clear that compositional variation has a much 
larger effect for the outer HZ limit, and for the higher mass stars in particular. 
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Habitable Zone Distance (AU) 

Fig. 9.— Inner and outer HZ edges (RGH and MaxGH) for all masses and compositions 
at TAMS. The elongated solid lines represent the 5 oxygen cases at Zq for each mass: 2 
enriched oxygen models (1.48 and 2.28 0/FeQ ), 2 depleted oxygen models (0.67 and 0.44 
O/Fee ), and 1 model representing the standard oxygen case (solar 0/Fe). The elongated 
dotted lines represent the end member values for the spread in oxygen abundance (0.44 and 
2.28 O/Fe© ) calculated at end member Z values (0.1 and 1.5 Zq ). From left to right are the 
highest to lowest opacity models. It is clear that compositional variation has a much larger 
effect for the outer HZ limit, and for the higher mass stars. Additionally, there exists a more 
exaggerated spreading trend for the lower mass stars at TAMS than occurred at ZAMS. 
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Age (Gyr) 



Age (Gyr) 


Fig. 10.— Inner (solid) and outer (dashed) edges of the HZ for three values of 0/Fe (at Zq ) 
for three masses: 0.5 Mq (top), 1 Mq (center), and 1.2 Mq (bottom). Each color represents 
a different O/Fe© value: black is solar 0/Fe, light gray is depleted (0.44 O/Fe© ), and dark 
gray is enriched (2.28 O/Fe 0 ). A 1 AU orbit is indicated by the dotted line. O abundance 
variations within a star can signihcantly affect MS lifetime as well as HZ distance. The inner 
radius is the Runaway Greenhouse case; the outer edge is the Maximum Greenhouse case. 
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Age (Gyr) 


Fig. 11.— Inner (solid) and outer (dashed) edges of the HZ at three Z values for three masses: 
0.5 Mq (top), 1 Mq (center), and 1.2 Mq (bottom). Each color represents a different Z value: 
black is Zq , light gray is depleted (0.1 Zq ), and dark gray is enriched (1.5 Zq ). A 1 AU 
orbit is indicated by the dotted line. Metallicity variations within a star can signihcantly 
affect MS lifetime as well as HZ distance. The largest effect is seen at 0.1 Zq . The inner 
radius is the Runaway Greenhouse case; the outer edge is the Maximum Greenhouse case. 
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Distance from Parent Star (AU) 


Fig. 12.— Inner (black) and outer (dark gray) boundaries of the HZ at the ZAMS (solid) 
and TAMS (dashed). The light gray shaded region is the CHZ, where an orbit would remain 
in the HZ for the star’s entire MS lifetime. This is for solar-composition stars for each mass 
in our range. For conservative HZ limits (RGH and MaxGH), low-mass stars have no GHZ. 
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Fig. 13.— Inner (black) and outer (dark gray) boundaries of the HZ at the ZAMS (solid) 
and TAMS (dashed). The inner edge 2 Gyr after the ZAMS and the outer edge 2 Gyr before 
the TAMS are indicated by dotted lines. The light gray shaded region is the GHZ 2 , in which 
an orbiting planet would remain in the HZ for at least 2 Gyr. This is for solar-composition 
stars for each mass in our range. The fraction of the total habitable orbits in the GHZ 2 is 
higher for the long-lived, low mass stars. 
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Fig. 14.— Inner (solid) and outer (dashed) edges of the HZ with cold starts prohibited 
for three values of 0/Fe (at Zq ) for three masses: 0.5 Mq (top), 1 Mq (center), and 1.2 
Mq (bottom). Black is solar 0/Fe, light gray is depleted (0.44 O/Fe© ), and dark gray is 
enriched (2.28 O/Fe© ). A 1 AU orbit is indicated by the dotted line. The inner radius is 
Runaway Greenhouse and the outer edge is the Maximum Greenhouse, at ZAMS value. 
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Distance from Parent Star (AU) 


Fig. 15.— Inner (black) and outer (dark gray) boundaries of the HZ at the ZAMS (solid) 
and TAMS (dashed). This is for solar-composition stars for each mass in our range. Here 
we equate the TAMS outer limit with the ZAMS outer limit in order to address the problem 
of cold starts. The thin-lined shaded region represents the range of orbits in which a planet 
would be in its host star’s HZ from the beginning of the MS, while the thick-lined shaded 
region represents the range of orbits that would be continuously habitable for the star’s 
entire MS lifetime (as in Figure [12]). If we attempt to understand a CHZ in this way, we see 
that the ideal habitable regions do not include the TAMS outer limit, where a planet could 
potentially start out very cold and only enter the HZ late in its lifetime. 
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Distance from Parent Star (AU) 


Fig. 16.— Inner (black) and outer (dark gray) boundaries of the HZ at the ZAMS (solid) 
and TAMS (dashed) for solar-composition stars for each mass in our range for the case. Here 
we represent the case in which cold starts are prohibited. The TAMS outer limit of the HZ 
is equated with the ZAMS outer limit. The shaded region now represents the range of orbits 
in which a planet would be continuously habitable for at least 2 Gyr and would also be in 
its host star’s HZ from the beginning of the MS. 
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